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Abstract. We compute the covariant three-point function near horizon-crossing for a sys- 
tem of slowly-rolling scalar fields during an inflationary epoch, allowing for an arbitrary 
field-space metric. We show explicitly how to compute its subsequent evolution using a co- 
variantized version of the separate universe or l SN } expansion, which must be augmented by 
terms measuring curvature of the field-space manifold, and give the nonlinear gauge trans- 
formation to the comoving curvature perturbation. Nonlinearities induced by the field-space 
curvature terms are a new and potentially significant source of non-Gaussianity. We show 
how inflationary models with non-minimal coupling to the spacetime Ricci scalar can be 
accommodated within this framework. This yields a simple toolkit allowing the bispectrum 
to be computed in models with non-negligible field-space curvature. 
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1 Introduction 



Considerable effort has recently been invested in the study of multiple-field models of infla- 
tion. There are three principal motivations: First, unless curvature couplings flatten their 
potentials at high energy density, the Standard Model has not produced scalar fields which 
can successfully inflate. This has stimulated the search for realizations of inflation in theories 
beyond the Standard Model. Second, single-field inflationary models often require super- 
Plan ckian field excursions. Unfortunately, well-rehearsed arguments suggest that we should 
not expect the scalar potential to be stable against renormalization-group running over such 
large distances in field space. Multiple-field models may evade this problem by allowing sub- 
Planckian excursions. Finally, interactions between several fields can give rise to observable 
non-Gaussianity. This may make some multiple-field models sufficiently predictive that they 
can be falsified by observation. 

It has been known for some time that, under slow-roll conditions, multiple- field models 
with canonical kinetic terms generate unobservable three- and four-point functions at horizon- 
crossing [1—4]. Observable effects can arise only later, from nonlinear processes operating on 
superhorizon scales. A reasonably clear picture has emerged in which these nonlinear effects 
can be understood as the deformation of a Gaussian probability distribution by the phase 
space flow associated with the theory [5-8]. 

Noncanonical kinetic terms offer new possibilities. In some theories the Lagrangian 
becomes an arbitrary function of the kinetic energy 2X = —d^ipd^ip. Where this yields a 
reduced sound speed for perturbations there can be a significant enhancement of the three- 
and four-point functions [9-15]. But this is not the only type of noncanonical Lagrangian. In 
many examples descending from our present ideas about physics at very high energies, includ- 
ing string theory and supergravity, the kinetic energy must be written 2X = —Gijdn<p d^ip , 
where Gjj is an arbitrary, symmetric function of the fields ip . The simplest example is the 
nonlinear er-model of Gell-Mann & Levy, originally introduced to describe spin-0 mesons. 

The matrix G/j can be thought of as a metric on the space of scalar fields and will 
generically exhibit nonzero curvature. Are there interesting enhancements of non-Gaussian 
effects in these models? Estimates of the three-point function have been made by a number 
of authors [14-18], but as yet no complete formalism exists which allows the bispectrum to 
be followed from horizon-crossing up to the time of observation. Moreover, Gong & Tanaka 
recently pointed out that the most widely-used formulation of nonlinear perturbation theory 
is not covariant [19]. They introduced a covariant description, to be discussed in §2.1, and 
constructed the action for fluctuations up to third-order in the scenario of Langlois et al. [14]. 
A similar argument was later made by Saffin [18]. Covariance is a convenience rather than 
a physical principle, so its absence does not invalidate earlier results. Nevertheless, it is a 
considerable convenience: there are subtleties associated with time evolution of the two-point 
function on curved field-space which are most clearly expressed in covariant language. These 
take the form of a contribution to the effective mass-matrix from the Riemann curvature 
tensor. In this paper we show that a similar phenomenon occurs for the three-point function. 

In flat field-space, time evolution of superhorizon modes may be taken into account using 
the 'separate universe' method [20-23]. This enables the time dependence of each fluctuation 
to be determined from the relative behaviour of separated spatial regions following slightly 
displaced phase-space trajectories. It can be effected using a Taylor expansion to compare 
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two solutions of the slow-roll equation. But in curved field-space we must be cautious when 
comparing the relative motion of neighbouring trajectories. In the analogous case of general 
relativity one would use the equation of geodesic deviation, or 'Jacobi equation'. In flat 
field-space this can be integrated to reproduce the Taylor expansion [8]. When promoted to 
curved field-space the Jacobi approach is automatically covariant and accounts naturally for 
time-dependent effects generated by the Riemann curvature, including its known contribution 
to the effective mass-matrix. Beyond linear order there are new contributions which influence 
the three-point function. These must appear in any correct formulation but are expressed 
most clearly and economically in terms of the field-space curvature. 

At present, the covariant approach cannot be used to generate observable predictions 
beyond the power spectrum. To do so would require a determination of the initial value of 
the three-point function near horizon-crossing, together with a prescription to evolve it into 
the primordial curvature perturbation. It is only the curvature perturbation which can be 
connected with observable quantities. Neither of these pre-requisites has yet been provided. 

In this paper we compute the initial value of the covariant three-point function at 
horizon exit and use the Jacobi approach to determine its time evolution. (Partial expressions 
for the noncovariant three-point function were given by Langlois et al. [14] and Renaux- 
Petel et al. [15].) As a concrete example we give the analysis for the a-model Lagrangian 
C = X + V, although our methods generalize to more complex cases. We show that this 
initial value can be smoothly connected to the subsequent Jacobi evolution. In particular, the 
evolutionary effects described above which depend on the Riemann tensor can be matched to 
new infrared divergences in the three-point function. We demonstrate this matching explicitly 
to subleading order in both time-dependent perturbation theory and the slow-roll expansion. 
A further benefit of the Jacobi approach is that time evolution can be computed very simply 
using ordinary differential equations. 

In §2 we specialize the results of Gong & Tanaka to the cr-model Lagrangian and obtain 
the action for fluctuations to third-order. In §3 we compute the corresponding two- and three- 
point functions near horizon-crossing in the spatially flat gauge. The two-point function has 
been known since the work of Sasaki & Stewart [24], but the computation of the covariant 
three-point function is new. In §4 we use the Jacobi method discussed above to compute the 
evolution of these correlation functions after horizon-crossing. 

In §5 we show that our results can be applied to models in which the scalar fields are 
coupled non-minimally to gravity by making a suitable conformal transformation [25]. Such 
couplings arise naturally in the low-energy limit of higher-dimensional theories including 
supergravity, string theory and Kaluza-Klein models [26-28], or as counterterms in curved 
spacetime [29, 30]. Finally, we conclude in §6. 

Notation. Throughout, we work in units where c = h = 1 and express the gravitational 
coupling in units of the reduced Planck mass, M~ 2 = 8ttG. Upper -case Latin indices /, J, 
K, . . . , label the species of scalar fields, and Greek letters label spacetime indices. We use 
a modified index convention for bilocal tensors ('bitensors'), to be described in §3.3. The 
covariant derivative compatible with the field-space metric G IJ is V/. For any tensor T... 
we write V/J 7 ... = J-.-.-j. Our sign convention for the curvature tensor is defined by the 
Ricci identity, [V/,Vj]Vr- = RijklV L - Finally, it is useful to define covariant versions of 
the derivatives with respect to coordinate time t, conformal time r\ = J dt/a and e- folds 
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N = jHdt as 

dt ' J ' 

2 The action and its perturbations 



^ = ^V„ 2>„ = ^V„ V N = %Vt. (1.1) 



Consider an inflationary epoch driven by Af scalar fields ip (with I = 1, 2, . . . ,Af), minimally 
coupled to gravity and self-interacting through a potential V(ip), 

S= l -Jd A x^ [M^R - GusTd^dy^ ~ W] ■ (2-1) 

As explained in §1, the field space is to be understood as a manifold with metric Gu(ip). 
This metric is used to raise and lower tangent-space indices, and distinguishes the model 
from canonical scenarios where Gu = Sjj. 

During inflation the ip 1 take approximately homogeneous but time-dependent values 
which we label (f) 1 (t). Small inhomogeneities around this homogeneous background are con- 
trolled by a perturbative expansion of (2.1). When quantized they will seed the primordial 
inflationary fluctuation. 

2.1 Covariant perturbations 

In what follows we perturbatively expand the action for these fluctuations to third-order. 
On curved field-space it is very helpful to arrange this expansion so that covariance is mani- 
fest. Each inhomogeneity is a coordinate displacement Sip 1 = y/(x, i) — ^(t), but for finite 
length this displacement does not lie in the tangent space at cf) 1 and therefore does not have 
a tensorial transformation law. To obtain a covariant description we must find an alterna- 
tive characterization which associates each displacement with a tangent-space vector. This 
construction was given by Gong & Tanaka [19] , whose method we briefly describe. 

We assume the field-space metric to be smoothly differentiable in the neighbourhood of 
the background trajectory. Within a normal neighbourhood, the points ip 1 (~x-,t) and (f) 1 (t) are 
linked by a unique geodesic, labelled by a parameter A. We adjust the normalization so that 
A = corresponds to the unperturbed coordinate c/) 1 and A = 1 corresponds to the perturbed 
coordinate tp 1 + Sip 1 . The coordinate displacement Sip 1 may be expressed as a formal Taylor 
series along this geodesic 



Sip 1 



dip 1 



dX 



+ ^7 



i dV 



a=o 2! dA^ 



1 d 3 ip 



a=o 3! dA3 



+ •••. (2.2) 

A=0 



Eq. (2.2) is independent of the normalization of A, so our particular choice was merely a 
convenience. The geodesic satisfies 

7^2 T d V , r l dip J dip K 

where T>\ = Q 7 V/ and Q 1 = d<^ / /dA|A=o- Using (2.3), the expansion (2.2) can be reorganized 
as a power series in Q 1 , yielding 

V = Q 1 - ^j K Q j Q k + 1 Km^jk - rj K , L ) Q J Q K Q L + ■■■, (2.4) 
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where the coefficients Tj K , Fj KL , . . . , are evaluated at A = 0. In flat field-space all terms 
but the first vanish and therefore Sip 1 = Q 1 . It was to achieve this correspondence that 
we adopted our normalization for A. Although (2.4) is not itself covariant, it can be used 
to exchange an expansion in Sip for an expansion in Q 1 . Since Q 1 does lie in the tangent 
space at (j) 1 , the perturbative expansion of any tensorial object will be manifestly covariant if 
expressed in powers of Q 1 . We label tangent-space indices at the perturbed position (j) 1 + Sip 1 
by primed indices I', J' and tangent-space indices at the original position <j) by unprimed 
indices. For any tensor Ti...j we obtain 

Fr...j. = Gp 1 ■ ■ ■ Gj, J f Jw| A=0 + ^/...j| A=0 + ^ V 2 F L ..j\ x=Q + "\ ( 2 - 5 ) 

where Gp 1 is the parallel propagator, which expresses parallel transport along the geodesic 
connecting (fr with cj) + Sip 1 . For details, see Poisson, Pound & Vega [31]. 

2.2 Gauge choice and infrared-safe observables 

Inhomogeneities in ip couple to gravity and therefore induce fluctuations in the metric. The 
description of this mixing is simplified using ADM variables [32] , in terms of which the metric 
can be written 

ds 2 = —N 2 dt 2 + hijidx* + N l dt)(dx j + N j dt), (2.6) 

where N is the lapse function and N l is the shift vector. Spatial indices are raised and 
lowered using the 3-metric hij, which has determinant h. In these variables, the action can 
be written 



s = i j d A x VK |m p 2 l iv,R( 3 ) + -^i /•:,,/•." - E 2 ) 



+ ^ 7 tt/ - Ndi^d 1 ^ - 2NV 



(2.7) 

where is the Ricci scalar built from the 3-metric and Eij is proportional to the extrinsic 
curvature of slices of constant t, 

2/w, = - N Aj - N j{i . (2.8) 

We have defined E to be the trace E l i, and a vertical bar denotes the covariant derivative 
compatible with hij. Finally, tt 1 = ip — N^ip y and an overdot denotes a derivative with 
respect to t. 

The lapse and shift appear in (2.7) without time derivatives. Therefore they are not 
propagating modes, but yield constraint equations. These determine N and N l as functions 
of the physical degrees of freedom, and also generate gauge symmetries associated with 
translations in time and space. After imposing the constraints and removing redundant 
gauge modes we are left with M scalar degrees of freedom plus two polarizations of the 
graviton which transform as a spin-2 mode, but we are free to choose how the scalar modes 
are divided between the ip 1 and hij by making gauge transformations. 

Infrared-safe observables. When selecting a gauge we must balance competing demands. 
First, consider a single-field model where N = 1. A common gauge choice is to arrange slices 
of constant t to coincide with slices of constant ip, leaving a scalar metric mode which we 
label C- This scalar mode is a perturbation to the volume element, £ = (1/6)5 In det h. The 
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disadvantage of this choice is that the calculation is long [33] . Multiple partial integrations 
are required to bring the action to a suitable form, simultaneously generating boundary terms 
which must be retained [34-36]. 

The advantage is a precise type of technical simplicity. With only a single field it is a 
classical theorem that £ is time independent [37, 38]. The same is true in the quantum theory, 
at least at tree-level, for the correlation functions of £ [39]. 1 If present, time dependence would 
manifest itself as a divergence in the vertex integrals appearing in each n-point function [43- 
45]. Each integral sums the amplitude for an interaction to occur at a specific time, and 
a divergence means only that interactions continue arbitrarily far into the future. At tree- 
level the great merit of £-gauge for single-field models is that the n-point functions are free 
of divergences: each vertex integral receives significant contributions only from interactions 
which occur near horizon-crossing. Physically, predictions for £ decouple from the infrared 
dynamics of the theory. 

For this reason we describe £ as 'infrared safe'. 2 In practice this means that subleading 
terms in the action map to subleading terms in each correlation function. This makes it 
simple to impose approximation schemes, such as the slow-roll expansion. 

The enormous convenience of infrared safety means that most single-field calculations 
have made use of this gauge. An alternative is to fix slices of constant t to carry a flat metric, 
leaving a perturbation in the scalar field value (p(x,t). In view of the discussion in §2.1 we 
denote the covariant representation of this perturbation Q 1 . Calculation of the lowest-order 
action for Q 1 is simpler than in the C-gauge [1]. But unlike £ the field perturbation is not 
infrared safe because the integrals which define its correlation functions receive contributions 
from all times, not just those near horizon crossing [43, 47-50]. Therefore they are sensitive 
to the infrared dynamics of the theory. As a result, subleading terms in the action can be 
enhanced by divergences and special action is required to deal with them. In a single-field 
model, Maldacena argued that they can be accounted for by a gauge transformation to £ just 
after horizon crossing [33]. 

Multiple- field models. In a multiple-field model the situation is more complex. Because £ 
is evolving, the divergences can no longer be accounted for using only a gauge transformation 
on a fixed time-slice. The solution proposed in Ref. [1] was to evaluate each n-point function 
in the spatially flat gauge a few e-folds after horizon crossing. This choice prevents enhanced 
subleading terms from spoiling the lowest-order slow-roll prediction, and therefore we can 
take advantage of the computational simplicity of this gauge. The price we pay is that some 
means must be found to express the correlation functions at a subsequent time in terms 
of their values near horizon crossing. However this is done, the result must equal what 
would have been obtained if we had been able to evaluate the original divergent integrals. 

: A similar statement can be made at one-loop level, ignoring internal graviton lines [40, 41]. Even if 
present, loop-generated time dependences are typically strongly suppressed [42]. Although they might be 
important to describe the evolution of correlations over very large time and distance scales, it seems probable 
they would be negligible for the description of a phase of observable slow-roll inflation. 

2 The terminology is borrowed from gauge theories, where an infrared-safe observable dominated by a hard 
subprocess occurring near energy E does not depend on the details of other processes (such as hadronization 
and confinement) occurring at energies much less than E. A similar discussion was given by Weinberg [46], 
who focused on the conditions under which (in our language) observables might be infrared-safe. (However, 
there is no reason of principle for any physical observable of interest to be infrared-safe.) 



- 6 - 



Therefore, in a perturbative expansion, it must reproduce the same divergences. For canonical 
fields on flat field-space, the nonlinear separate universe or formalism introduced by 
Lyth & Rodriguez can be used for this purpose [23]. 

On curved field-space we expect new divergences involving the field-space Riemann cur- 
vature [19] , and we must be aware that these will modify the time-evolution of the correlation 
functions. We will argue in §4 that these can be understood as an analogue of the geodesic- 
deviation effect for freely-falling observers in curved spacetime, and show how they can be 
incorporated in a new version of the l SN } formula. 

The conclusion of this discussion is that we are still entitled to choose the spatially 
flat gauge in order to simplify the calculation. However, we must take care to study the 
effect of divergent terms. Although we will carry the QFT calculation of the three-point 
function only as far as horizon-crossing, where the divergences are harmless, it is important to 
check carefully that whatever technique we employ to account for time dependence correctly 
reproduces all these divergent pieces. 



3 Two- and three-point functions of the fluctuations 

Flat gauge calculation. With this in mind, we specialize to the spatially flat gauge where 
hij = a 2 5ij. To expand (2.1) to third-order in Q 1 we need only solve the constraints for the 
first-order components of N and N l [11, 33]. The shift vector can be decomposed into 
irrotational and solenoidal parts, yielding N l = <9*# + where <9j/3' = 0. The first-order 
solenoidal component (3\ appears uncoupled to the other perturbations in the second-order 
action S^), yielding the constraint equation (3\ = 0. The remaining metric perturbations can 
be expanded in powers of Q, 



N = 1 + ai + a 2 + 

= #1 + #2 + ■ • • • 



(3.1) 



In what follows we determine these metric fields and use them to obtain the two- and three- 
point functions for the Q 1 . 

3.1 Linear order 

At linear order we find 



5(i) 



<Tx I a 



2 9IC 



V 



V t {a 3 </>/) + a 3 Vr 



Q 1 



(3.2) 



where we have integrated by parts and removed total derivatives. The background field 
equations follow after varying this action with respect to a\ and Q 1 , 



3M 2 ,H 2 = l -4>^ + V, 



V t <j>! + = -Vt. 



(3.3) 
(3.4) 



The slow-roll regime in curved field-space was discussed by Sasaki & Stewart [24] and later 
Nakamura & Stewart [49]. Inflation occurs when e = —H/H 2 < 1. For a successful phe- 
nomenology we also require that inflation is sufficiently prolonged. Therefore e should not 
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change significantly in a Hubble time. To satisfy these requirements we must choose 

1 ¥h , _ 1 dine 2 ^V^j 

e= 2M% l lP <<1 and ^H^T = H^ + 2e<<1 - (3 ' 5) 

The condition i] <C 1 requires the tangential component of the acceleration vector H -1 !)^ 1 
(measured in Hubble units) to be much smaller than the tangent vector to the trajectory. 
One can verify that the slow-roll equation 

3#0 7 + Vj fa (3.6) 

gives a self-consistent realization of these conditions if the potential is sufficiently flat. De- 
tailed conditions are given in Refs. [24, 49]. However, global flatness of the potential is not 
necessary. This possibility has received recent attention [51-54]. In this paper we do not 
assume any particular properties of the background theory, except that it realizes an era of 
inflation with e <C 1 and in which H varies smoothly during horizon exit. Such scenarios 
are not the only cases of interest, but a dedicated analysis is required where the background 
evolution exhibits a feature [34, 55]. 

When we compute the two- and three-point functions we will do so only for field-space 
directions which are light during horizon exit. Therefore our results will not apply to any 
heavy directions generated by a steep potential orthogonal to the inflationary trajectory. 
(However, when estimating the magnitude of terms in the second- and third-order actions, 
we quote powers of ft /H to emphasize that individual components of this vector are not nec- 
essarily of order e 1//2 .) This is sufficient to estimate the statistics of the primordial curvature 
perturbation in an approximation where the fluctuations in massive directions at horizon 
exit are neglected. In simple models this is acceptable because large masses rapidly drive 
any fluctuations to extinction. In more complex models it has been suggested that modest 
corrections can occur where the phase-space flow drives power from massive modes into the 
curvature perturbation before decay [56-58]. To capture these effects would require an ex- 
tension of the formalism of §§3.2-3.3 used to compute initial conditions, although we expect 
that the subsequent transfer of power would be correctly described by the method of §4. 

3.2 Second order 

Expanding the action to second order, performing multiple partial integrations and removing 
total derivatives, we find 



S {2) = ^J d 4 xa 3 {ai - 6M^H 2 ai + fatfctx - 2<t> I V t Q I - IV jQ 1 



(3.7) 



+ Rkijl<P K <P L Q I Q J + VtQjVtQ 1 - h^diQfdjQ 1 - ^/jQ 7 Q J }. 



The momentum and energy constraints can be obtained by varying the action with respect 
to ot\ and We find 

2M^Ha x = faQ 1 , (3.8) 

-2M p 2 ^5 2 i?i = 6MpiH 2 a\ - a^tf + faVtQ 1 + VjQ 1 . (3.9) 
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In writing Eqs. (3.7) and (3.9) we have used the spatial Laplacian d 2 = <9j<9j, with the 
convention that spatial indices which are both written in the covariant position are summed 
using the flat Euclidean metric 5 13 . Employing Eqs. (3.8) and (3.9) we can eliminate the 
metric perturbations in 5(2) to obtain 

S(2) = \J ^xa^VtQiVtQ 1 - Wd&idjQ 1 -MijQ t Q J }, (3.10) 
where the symmetric mass matrix «Mjj satisfies 

Mu = Vu - R L UM<p L <j> M - j^V t (jjHh^J ■ (3.H) 

This is identical to the canonical case except for the covariant derivatives and the term 
involving the Riemann tensor, which was first obtained by Sasaki & Stewart [24]. (See also 
Nakamura & Stewart [49] and Gong & Stewart [59, 60].) We will find similar terms in the 
third-order action (3.14) below. Their meaning is not immediately clear because they imply 
that promotion to curved field-space requires more than 'minimal coupling' to curvature. In 
the mass matrix, the Riemann term changes the effective mass of modes orthogonal to (j) but 
also alters the coupling of these modes to each other. We discuss these effects in more detail 
in §4. 

Power spectrum at horizon-crossing. To compute the power spectrum and all higher 
n-point functions we must use the 'in-in' formulation of quantum field theory, which entails 
doubling all field degrees of freedom. For details, we refer to the literature [4, 46, 61, 62]. 
In curved field-space an extra complication is caused by the necessity to give each n-point 
function the correct tensorial transformation properties. As we describe in §3.3 below, this 
is enforced by transport of the tangent-space basis along the inflationary trajectory. 

The power spectrum at lowest order in M.u was calculated by Sasaki & Stewart [24] 
and can be obtained from (A. 10) or (A. 11). Taking the equal-time limit in either equation, 
it follows that the two-point function evaluated a little later than horizon-exit satisfies 

(Q^OQ'Cka)) = (2^) 3 ,5(k 1 + k 2 )^G /J . (3.12) 

This estimate becomes valid when the decaying power-law terms in (A. 10) and (A.ll) have 
become negligible [63]. However, because Q is not infrared safe there are growing terms 
at subleading order [48, 49]. Therefore (3.12) soon becomes untrustworthy, and its validity 
extends only for a very narrow range of e- folds. This is the problem of enhanced subleading 
terms discussed in §2.2. In Eq. (4.16) we use the separate universe method to give a more 
accurate expression which remains valid until late times. 
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3.3 Third order 



The third order action is 

S {3) = - f d A xa 3 \6M^H 2 al + iM^afd 2 ^ -^^-(d i d j 'd 1 did j 'd l - d^d 2 ^ 
2 J I a a V 

- fatfa\ + 2o? 1 j> I V t Q I + -^a^d^diQ 1 - a\R L ^ IJ ^ M 4>l4>mQiQj 

or 

- aJVtQ^tQ 1 + \diQtdiQ 1 ) - ^d^VtQj^Q 1 + ^Ri i j K )J L ^tQ I Q J Q K 

\ ex / a o 

+ \R{i\lm\J;K)^ L 4> M Q I Q J Q K ~ lv. (IJK) Q I Q J Q K - %j)aiQV}- 

(3.13) 

We have indicated the symmetric index combinations picked out by each product of the 
Q 1 . 3 The lapse and shift can be eliminated using (3.8)-(3.9). The resulting expression is 
exact and does not invoke an expansion in powers of slow-roll. Nevertheless, the argument 
of §2.2 implies that we need only compute lowest-order contributions provided we evaluate 
the three-point function near horizon-crossing. In this region infrared divergences cannot 
enhance sub leading terms. Focusing only on the low-order contributions and rewriting in 
terms of conformal time, we find [19] 



S( 3 ) 5 



,2 



+ WJ72-E^ d ^ 2 ' D vQi d iQj' D vQ J ( 3 - 14 ) 

2ivt pl n 

+ ^Ri(jK)J L v ri Q I Q J Q K + ^R { i\lm\j;k^ l ^ m Q i Q j Q k Y 

The first two lines of (3.14) are a covariantization of the action obtained in flat field-space [1]. 
We will sometimes describe them as the 'canonical' terms. Subleading corrections begin at or- 
der 0((j)/H) 3 and are negligible unless enhanced by divergences. The third line includes terms 
containing the Riemann tensor, analogous to the Riemann term in the mass matrix (3.11). 
These were first obtained by Gong & Tanaka [19] and are a new feature associated with the 
curvature of field space. Unlike the lowest-order 'canonical' terms they produce divergences. 
To track the influence of these as clearly as possible we have retained Riemann terms up to 

Although operators involving the curvature at 0((p/H) 2 are subleading, we can reliably 
compute contributions to the three-point function at this level because Eq. (3.11) shows that 
next-order corrections from the propagator are themselves 0(<p/H) 2 , and therefore enter the 
three-point function only at 0(4>/H) 3 . The same is true for corrections from the scale factor 
and Hubble rate. After expanding around the horizon-crossing time for a fiducial scale k^, 
the explicit 0(<j)/H) 2 term in (3.14) is accompanied by one further contribution at the same 
order from the time-dependence of Ri[jk)l- [See Eqs. (A. 15) and (A. 22).] We will retain 
both these terms when estimating the three-point function. The advantage of doing so is 



3 Our symmetrization conventions are 2A^ IJ - j = Azj + Aji and GA(i JK ) = Auk + {5 perms}. Bars delimit 
indices excluded from symmetrization. 
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that we can perform a more stringent test of our matching to the superhorizon evolution 
in §4. 

Three-point function. In Appendix A we calculate the contribution to the three-point 
function from each of these operators. As explained in §3.1, our computation applies only to 
light field-space directions for which Mu can be neglected. 

To express the result it will be necessary to perform parallel transport along the infla- 
tionary trajectory. This can be accomplished using the parallel propagator, 

M' 



U T i = Pexp j\ v T r M , N/ ^— Wi, (3.15) 



where the integral is computed along the inflationary trajectory between conformal times 
a and r. The symbol 'P' denotes path ordering along this trajectory. (More details and 
sample calculations are given in Appendix A.) We refer to H 1 i as the trajectory propagator. 
In writing (3.15) we have adopted a notation in which the index / is associated with a basis 
for the tangent space at time r, whereas the index % is associated with an independent basis 
for the tangent space at a. Primed indices label the tangent space at a time corresponding 
to the integration variable rj . The rightmost N' index lies at time a and contracts with 
G n i. Therefore Il^j is a bitensor: it is an object with mixed indices [31], transforming like 
a contravariant vector at r and a covariant vector at o. In what follows we set a to be the 
horizon-crossing time for the reference scale A;*. 

When computing n-point functions we typically measure time in e-folds of expansion, 
evaluating each n-point function at iV = — ln|fe*r| e-folds since the fiducial scale k* passed 
outside the horizon. With these conventions we find 

(Q I (k 1 )Q J (k 2 )Q K (k 3 )) = (2vr) 3 ,5(k 1 + k 2 + ^J^Il 1 \U J jU K k A^ k {N). (3.16) 

The bitensor A^ k (N) transforms as a scalar at time N and a rank-three tensor at the time of 
horizon-crossing for fc*. However, the three-point function (Q 1 Q J Q ) transforms as a rank- 
three tensor at N. To match these transformation properties, the propagators IT^, II J j and 
H k are necessary. In Appendix A we explain how they arise in terms of Feynman diagrams. 
We set kt = k\ + &2 + &3 to be the total scalar momentum and define k 2 = Ylii<j hkj. Under 



4 In (2.5), to match the notation of Ref. [31], we denoted the parallel propagator evaluated along a geodesic 
as Gii 1 . We reserve the symbol II to indicate parallel propagation along an inflationary trajectory. 
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the same conditions which were required for validity of (3.12), we find 



A ijk (N) 



M 2 H. 



2 MM + h k2 . k3 



kt 7E 



+ 3 R * Hi 



1 R (i\mn\j;k) <j>m <j>n 
3 H* 

_ 4 R i(jk)m;n Ym Tn 

3 H* Hi. 



In 



+ ^3 



7E 



Mr 



2 V 



7T 



7E+ 12 



2 7E + In 2- 



ln^ 



fci fc t In — + 7e - 1J — (je + In 



+ cyclic, 



(3.17) 



where 7e ~ 0.577 is the Euler-Mascheroni constant. The first line of (3.17) covariantizes 
Eq. (69) of Ref. [1], and subsequent lines arise from the Riemann terms in (3.14). The sum 
over cyclic permutations includes the two permutations generated by simultaneous exchange 
of {i, ki} with either {j, k.2} or {k, h^}. 



4 Evolution after horizon exit 



Terms in (3.17) involving N are divergent in the late-time limit r — > and are responsible 
for spoiling infrared safety, as described in §2.2. They generate time evolution after horizon 
exit [43-45] and rapidly invalidate the expressions derived in §3. In this section we show 
that the divergent Riemann curvature terms in (3.17) have a geometrical origin and explain 
how their contribution can be taken into account using a covariant version of the 'separate 
universe' method [8, 23]. The terms diverging like a single power of N are 

a?"Jl = ^Nkf ( R(i\™\m^k _ + cyclic . (41) 

3 V H H H J * 

Eq. (3.17) also contains an explicit double logarithm (a term proportional to iV 2 ) 

4?og = ~\ N2k l (V^"^) + cyclic. (4.2) 

In addition, it implicitly contains terms of all powers in N from higher-order corrections we 
have not evaluated. 

Separate universe approach. Evolutionary effects on superhorizon scales can be under- 
stood using causality and classicality. After smoothing over small-scale structure, widely 
separated regions locally evolve like an unperturbed or 'separate' Friedmann universe [20- 
22]. Lyth & Rodriguez [23] extended this method to n-point functions for n ^ 3. In their 
formulation, the principal tool was a Taylor expansion of the background solutions in small 
deviations from a chosen initial condition. 
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In Ref. [8], the separate universe approach was applied without directly invoking this 
Taylor expansion. The inflationary trajectories in phase space are interpreted as the integral 
curves of a flow, and can converge or disperse. In Ref. [8] it was shown that the Taylor 
coefficients used by Lyth & Rodriguez can be understood as (derivatives of) the Jacobi fields 
which describe this dispersion. The growth and decay of fluctuations, and the processes by 
which power is transferred between modes, can be understood as the local dilation, shear 
and twist of a narrow bundle of trajectories. 

4.1 Jacobi equation for separate universes 

The Jacobi method provides a simple way to implement the separate universe approach in 
curved field-space. Consider two separate universes, with slightly displaced initial conditions, 
which correspond to neighbouring trajectories on phase space. The displacement between 
these trajectories can be described covariantly using a tangent-space vector Q (the 'con- 
necting vector') as described in §2.1. 

Deviation equation. Each universe evolves according to the field equation 

IpV+(l-f)zW = « J . (4-3) 

where uj = —Vj/3H 2 . Under the slow-roll approximation the acceleration term 'Dj^ip 1 is 
negligible along each trajectory. In flat field-space this means that the change in acceleration 
term between neighbouring trajectories also contributes at higher-order in slow-roll. On 
curved field-space this is no longer true because derivatives do not commute. Therefore we 
must retain the acceleration term when studying how trajectories disperse. 

The evolution of Q 1 can be determined by making a Taylor expansion of (4.3) along 
a geodesic connecting the adjacent trajectories, as in §2.1. To describe evolution of the 
two- and three-point functions we require this expansion up to second-order. Dropping the 
explicit 0(e) term, which can contribute only at higher order in the slow-roll expansion, and 
discarding a common factor of the parallel propagator we find 

V x + \v\ \ (\vW + T> N <A = u 1 ,,^ + l -u I , JK Q J Q K . (4.4) 

Using the Ricci identity to commute T>\ with T>n and employing the Bianchi identities to 
symmetrize resulting curvature terms, we conclude that Q 1 evolves according to the Jacobi 
or 'deviation' equation 

VnQ 1 = w'jQ- 7 + \w\ JK) Q J Q K + • • • , (4.5) 

where the coefficients w 1 j and w 1 \jk) satisfy 

1 cj) L <j) M , . 

WIJ = «(J;J) + 3 RL (IJ)M-y--fti ( 4 - 6 ) 

1/ $ L $ M </> L \ , 

Wl(JK) = U(I;JK) + ^[ R (I\LM\J;K)Jj^- ~ 4Hl(JK)LJj J ■ I 4 -?) 
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To obtain (4.6)— (4.7) we have imposed the slow-roll approximation to determine T^Q 1 in 
terms of V^Q 1 . As usual, the background trajectory is denoted ^(t). All curvature quanti- 
ties and derivatives of u 1 are evaluated on this trajectory and therefore powers of slow-roll can 
be counted in the usual way. Because we have used the slow-roll approximation, Eqs. (4.6)- 
(4.7) are trustworthy only to lowest slow-roll order in the derivatives of u 1 , and to 0(<p/H) 2 
in terms multiplying the Riemann tensor and its derivatives. This accuracy is sufficient to 
make a comparison with the divergent terms retained in §3.3. 

Although both terms in (4.6) are automatically symmetric under exchange of IJ we 
have indicated this explicitly. However, wijk is symmetric only on JK. This is different to 
the case of flat field-space, where terms involving the Riemann tensor are absent and each 
coefficient on the right-hand side of the Jacobi equation is always a symmetric combination 
of partial derivatives. When writing w 1 ijk) we a dd redundant brackets to emphasize this 
symmetry. 

Time-evolution operators. The Jacobi equation (4.5) is a first-order differential equa- 
tion, and therefore its solution can be expanded in powers of the initial conditions Q™, 5 

Q 1 = T'mQT + \T\ mn) QTQ T : + ■■■■ (4.8) 

To write (4.8) we have used the index convention introduced in §3.3. The fluctuation Q 1 is 
evaluated at some late time A, and its index I transforms as a contravariant vector in the 
tangent space at this time. Conversely, Q™ and its index m transform as a contravariant 
vector at an earlier time A*. Like the trajectory propagator (3.15), the coefficients T 1 m and 
T ( mn ) are bitensors. In particular, T m transforms like a contravariant vector on I and a 
covariant vector on m, whereas T 1 i mn \ transforms like a contravariant vector in the tangent 
space at A and a covariant rank- two tensor in the tangent space at A* . The initial conditions 
require T 1 m = 5^ and T 7 ( mn ) = when N = N*. 

Eq. (4.8) is a solution to the Jacobi equation (4.5) provided the T coefficients satisfy 

VnT^ = w I jT J m , (4.9) 
VnT 1 ^ = «/jT J (mn) + w\ JK) T J m T K n . (4.10) 

(Recall that A in the derivative T>n is not a field-space index, but the number of e-folds.) 
Eq. (4.8) can be summarized as the Taylor expansion of Q 1 in terms of its value at some 
earlier time A*, and defines a 'separate universe' approximation for curved field-space. We 
describe the coefficients T 1 ( m ... n \ collectively as 'time-evolution operators'. They are covariant 
analogues of the coefficients dcf) 1 /d(j)™ (and its higher derivatives) which occur when applying 
the separate-universe method in flat field-space [65, 66]. We could obtain them by solving 
ip 1 and using (2.5) to compute its derivatives with respect to the initial conditions, 6 but in 
practice it is much easier to integrate (4.9)-(4.10) directly. 

Divergences. We now show that (4.5)— (4.7) reproduce the divergences of (3.17). The 
argument is similar to that of Zaldarriaga [43]. Solving Eqs. (4.9)-(4.10) perturbatively 

5 The quantities T 1 m and T 1 mn were written F 1 m and F 1 mn in Refs. [8, 64]. In this paper we reserve T to 
mean the Levi-Civita connection compatible with the field-space metric G IJ . 

6 To reproduce all information in (4.6)-(4.7) it would be necessary to retain the connection term from V^tp 1 . 
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yields a power series in N. The lowest-order terms are 



T 1 —Tl 1 4- TT 7 ■ 

- 1 m — LL m \ L1 - i 



W 



N + -T1 1 



+ V N w l 



N 2 + 



(mn) 



u 1 



IV (mn) 



N 



(4.11) 
(4.12) 



+ 1* 



^Nw\ mn ) + W l k W k {ynn ) + w\ mk) W k n + w\ nk) W k 7 



N 2 + 



where N = — In |fe*r| as above. 

Eq. (4.11) shows that the time-evolution operator T 1 m can be understood as a mod- 
ification of the trajectory propagator to include the effect of time-dependence along the 
inflationary trajectory in addition to parallel transport. This follows because the trajectory 
propagator II 7 m satisfies (4.9) with w 1 j = 0. 

At linear order in N, the two- and three-point functions following from (4.11)-(4.12) 

are 



NHi 



(Q 7 (ki)Q J (k 2 )) 5 (2n) 3 5(k 1 + k 2 )-^Il I l Il J j w i J, 



(4.13) 



NHt 



(Q / (k 1 )Q J (k 2 )Q^(k 3 )) D (27Tf5(k 1 + k 2 + k 3 )—-^U^U J ,U 



1 ^j^kw^kf + cyclic. (4.14) 



The symmetry properties of w !J and w l ^ k ^ ensure that these expressions are symmetric under 
simultaneous permutations of the indices I, J, K and their associated momenta ki, k 2 , k 3 . 
At quadratic order in N and lowest order in slow-roll there is a contribution only from the 
T>j^w l i mn ) term in (4.12). This gives 



(Q I (k 1 )Q J (k 2 )Q K (kg)) D (2vr) 3 ( 5(k 1 + k 2 + k 3 ) 



N 2 Ht 
4IM 



x n^n^n 



J TfK. 



± R mm;n^<tA fc 3 + cydic 



(4.15) 



It can be checked that (4.13) reproduces the divergence in the two-point function (in- 
cluding the term involving the Riemann tensor) found by Nakamura & Stewart [49]. Com- 
paring (4.14) and (4.7), it can also be checked that the terms in w % ^ k ^ involving the Riemann 
tensor reproduce the divergence in (4.1). It was to enable a nontrivial check of this matching 
that we elected to keep divergences up to 0(4>/H) 2 in the Riemann-tensor terms of (3.14). 
Finally, comparing (4.15) and (4.2) it can be checked that the lowest-order double-logarithmic 
divergence is also correctly reproduced. At the accuracy of our present calculation it is not 
possible to check whether the divergences proportional to uu.jf.) also agree. These are a 
covariantized version of the divergences which appear in flat field-space. Since the same is 
true for the 'canonical' operators in the Lagrangian (3.13) we should expect agreement. 

In principle, higher-order terms in N and (f>/H could be retained in the perturbative 
expansions (4.11)-(4.12), which would enable a check of matching at all orders. Although such 
a check would be interesting, the matching at single- and double-logarithm order provides 
no reason to believe it would fail. We will not attempt it in this paper. 



7 In Refs. [49, 59, 60] the factors of H 1 i were omitted. 
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4.2 Transport equations 

The time-evolution operators enable us to determine each n-point function after horizon exit. 
Translating the formulae of Lyth & Rodriguez, we obtain 

(Q / (k 1 )Q J (k 2 )) = T I m T J n (Q m (k 1 )Q n (k 2 )), (4.16) 

and 

(Q / (k 1 )Q J (k 2 )Q A '(k 3 )) = T / i T J m T A ' n (g i (k 1 )Q m (k 2 )Q"(k 3 )), 

+ T\ lm) T K r T J s J 7S3 (Q'(ki - q)Q r (k 2 )),(Q m (q)g s (k 3 )), + cyclic, ^ 

where 'cyclic' denotes the two permutations of the second line in (4.17) obtained by exchang- 
ing {/, ki} with {J, k 2 } or {K, k 3 }. When there is no time evolution and w 1 j = 0, Eq. (4.17) 
reproduces (3.16). In combination, Eqs. (3.17), (4.6)-(4.7), (4.9)-(4.10) and (4.16)-(4.17) 
constitute the principal results of this paper. For comparison with observation it only re- 
mains to make a gauge transformation from Q 1 to the curvature perturbation £, for which 
we supply the relevant formulae in §4.3. 

Up to this point we have worked in a holonomic frame derived from the field-space 
coordinates, but other possibilities exist. Since an n-point function of the Q m evaluated 
at iV* transforms as a rank-n tensor in the tangent-space at time iV*, Eqs. (4.16)-(4.17) 
are manifestly covariant. As a result, we are free to select a basis for the tangent space 
independently at the early and late times iV* and N. 

Equations for two- and three-point functions. The approach given above is simple 
and emphasizes its similarity with familiar 6N methods, but it is also possible to write 
evolution equations for the n-point functions directly. In Ref. [8] these were described as 
transport equations. 

We write the two-point function as 

(Q 7 (k!)Q J (k 2 )) = (2^) 3 5(kx + k 2 )|^, (4.18) 

where Y* IJ is symmetric. The amplitude of the local mode of the three-point function can be 
parametrized 



(Q I (k 1 )Q J (k 2 )Q K (k 3 )) = (2vr) 3 ,5(k 1 + k 2 + k 3 ) 



a I(JK) a J(IK) a K(IJ) 



(4.19) 



2 3 1 3 1 2 

Direct differentiation followed by use of (4.5) yields the equations 

V N E IJ = w/ L £ iJ + w J L Z LI + • • • , (4.20) 
V N a I(JK) = w T L a L W + w\a 1( - LK ^ + w K L a^ + w\ LM) Y> LJ Y> MK + • • • , (4.21) 

where the omitted terms involve higher-order correlation functions and are negligible in 
typical inflationary theories. Following the method described in Ref. [8] it can be verified 
that (4.20)-(4.21) reproduce Eqs. (4.16)-(4.17). 
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4.2.1 Interpretation of Riemann terms 

It is now possible to understand the significance of the interactions in (3.17) mediated by the 
Riemann curvature, and the infrared divergences to which they give rise. They correspond 
directly to those terms in the separate-universe Jacobi equation (4.5) which measure deviation 
between nearby trajectories. The derivation of §4.1 makes clear that this effect is entirely 
analogous to geodesic deviation between freely-falling observers in curved spacetime. 

These new sources of time dependence arise mathematically from tidal effects in field 
space. Their physical meaning can be understood as follows. An initial perturbation 
generically represents a mix of adiabatic and isocurvature fluctuations. The isocurvature 
fluctuations differentiate between 'separate universes', and correspond to a choice of infla- 
tionary trajectory measured from the fiducial trajectory at Q 1 = 0. As these trajectories flow 
over field-space the proper distance between them will vary depending on the metric Gij. 
Each trajectory evolves independently in the absence of gradient terms which couple spatially 
separated regions. Therefore the isocurvature component of Q 1 represents a displacement to 
a fixed trajectory and must respond to this variation. 

These effects generate new sources of dilation and shear which contribute to the redis- 
tribution of power between adjacent trajectories. In addition, the nonlinear terms in w (jk) 
generate new sources of non-Gaussianity. The precise form of these nonlinear terms is unfa- 
miliar only because they are irrelevant in familiar applications of geodesic deviation — such 
as the focusing and defocusing of a bundle of light rays, where the connecting vector can be 
taken to be infinitesimal. 

The Riemann tensor is antisymmetric on its first and second pairs of indices. Since <j> is 
proportional to the tangent to the curve, we conclude that the Riemann contribution to w j 
is zero when either index is aligned with the adiabatic direction. Gong & Tanaka emphasized 
that this leads only to new couplings between isocurvature modes [19]. Eq. (4.20) shows that 
these couplings influence how the isocurvature modes share power between themselves, but 
do not cause power to flow between the isocurvature and adiabatic directions. Such a flow 
must be mediated by the potential through un. jy 

In the special case where the trajectory is an exact geodesic its tangent vector is parallel- 
transported proportional to itself. In this case, the adiabatic mode decouples completely and 
no power flows to or from it. 

4.3 Gauge transformation to the curvature perturbation 

For comparison with microwave background experiments or galaxy surveys we must compute 
the n-point functions of the primordial curvature perturbation, £. This is a measure of the 
local excess expansion between a spatially flat hypersurface and a uniform density hyper- 
surface with which it coincides on average. In curved field-space the computation can be 
performed economically using the method of Ref. [64]. We expand N as a function of the 
density p. Taking Ap to be the displacement from a point of fixed density p c to an arbitrary 
initial location, we find 

To determine the variation of (4.22) under a change in the initial location we expand along 
a geodesic, as in §2.1, along which both Ap and the differential coefficients will vary. The 
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variation of Ap satisfies 



ff(Ap) = -V-iQ 1 - \v.,uQ I Q J + ■■■ . (4.23) 



Therefore, up to second order, we can express £ as 



C = S(AN) = NjQ 1 + ^N IJ Q I Q J + • • • . (4.24) 



The coefficients Nj and Njj satisfy 



where 



dN 

Nj = ~V ]It (4.25) 
dp 

= ~V iZJ + ^VjVj + ^(Aj.j + Aj.j), (4.26) 



Al = ^-(y^ r ' Kv " K ' (427) 



dN 1 V 



dp ' M^V^Vj 

<^1 = _2 L_ , J v _ 

dp 2 M?j V'^V-i M^ y (V'JVi) 3 ' JK ' 



(4.28) 
(4.29) 



V y,J 

<5iV coefficients. Eqs. (4.25)^(4.26) are defined at a single point in field space; they are 
not bilocal in the sense of the coefficients T m and T ( mn ). We can obtain analogues of these 
bilocal coefficients using (4.8) to relate the Q 1 to their values at horizon-crossing. This yields 

C(N) = N m Q™ + ^N imn) QTQ: + • • • , (4.30) 

where N m and iV7 mn \ transform as scalars in the tangent space at N, and (respectively) rank- 
one and rank-two tensors in the tangent space at JV*. We describe them as l 5N coefficients'. 
They satisfy 

N m = NjT 1 ^ (4.31) 
N (mn) = N!T\ mn) + NjNjT^T^. (4.32) 

It follows that the two- and three-point functions of £ are determined by 

(C(ki)C(k 2 )) = 7V m iV n (Q m (k 1 )Q n (k 2 ))*, (4.33) 

and 

(C(k 1 )C(k 2 )C(k 3 )) = N l N m N n (Q l (k 1 )Q m (k 2 )Q n (k 3 )) Sf 

+ N {lm) N r N s J (Q'(ki - q)Q r (k2)>*<Q m (q)Q s (k 3 ))* + cyclic, ^ 

where 'cyclic' indicates the usual combination of permutations, as in Eq. (4.17). 

Eq. (4.32) implies that N^ mn ^ is the covariant derivative of N m with respect to a change 
in the initial conditions. Therefore (4.30) agrees with the covariant 5N expansion discussed 
by Saffin [18]. A similar expansion has already been used by Peterson & Tegmark [17]. 
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Observable quantities. The statistical properties of £ are typically expressed in terms of 
its spectrum and bispectrum 



(C kl Ck 2 ) = (27r) 3 <5(k 1 + k 2 )P f (fci), (4.35) 
(C kl Ck 2 Ck 3 ) = (27r) 3 (5(k 1 + k 2 + k 3 )B c (h,k 2 , h). (4.36) 

Constraints on the spectrum are expressed in terms of a dimensionless quantity Vq 

k 3 H(k) 2 
r C (k) = T^PdV = N m N n G mn (k)^-L, (4.37) 

where the argument k denotes evaluation at the horizon-crossing time for k. 

The bispectrum contains considerably more information. Its terms can be divided into 
two types: first, those which come from interference effects involving decaying modes near 
horizon exit; and second, those arising from interactions between only growing modes far 
outside the horizon. The first type can have arbitrary dependence on the A;- modes k±, k 2 , 
&3 which appear in the three-point correlation function. The second type appear only in 
the 'local' combination (kfk 3 ,)" 1 or its permutations. With canonical kinetic terms, applying 
global slow-roll conditions to the potential and assuming that only light fields contribute to 
C, Lyth & Zaballa [67] showed that the result of Ref. [1] implies only the 'local' bispectrum 
shape can be observable. 

With a nontrivial field-space metric, interactions involving the Riemann tensor may 
change this conclusion. Inspection of (3.17) shows that these interactions contribute terms 
of both types. Depending on the field-space curvature it is possible that 'nonlocal' contribu- 
tions in (3.17) could be enhanced, but to determine whether this happens would require an 
extension of the analysis in Refs. [67, 68]. We leave this interesting question for future work. 
On the other hand these terms certainly modify the evolution of the amplitude of each local 
shape, as discussed in §§4.1-4.2. 

Where the Riemann curvature is sufficiently small to make nonlocal contributions neg- 
ligible, Eq. (4.34) yields an analogue of the familiar '<W formula for the amplitude of the 
local bispectrum, 

f locals 5 N m N n N(™) 

/nl ~6 (N r Nr)2 ' ^.38) 

where N m and Nt mn \ were defined in (4.31)-(4.32). In this case, Eq. (3.17) can be interpreted 
to mean (as in the case of canonical kinetic terms) that the bispectrum generated at horizon 
exit is negligible: subsequent time evolution is necessary to generate an observable non- 
Gaussian signal. 

4.3.1 Backwards formalism 

To track the evolution of mixed two- and three-point functions for the complete set of fluctu- 
ations, including isocurvature modes, it is necessary to solve for all components of T 1 ' m and 
I ( mn ). In an A/"-field model, the linear coefficient T 1 m has A^ 2 independent components and 
the quadratic coefficient T 7 ( mn ) has Af 2 (Af + l)/2 independent components. But for some 
purposes we may require only the two- and three-point functions for £ given by (4.33)-(4.34). 
If so, we can reduce the computational complexity by tracking only the N components of N n 
and the M{M + l)/2 components of N^ mn y 
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An evolution equation for N n was given by Yokoyama, Suyama & Tanaka [69, 70]. As- 
suming the local shape dominates, only the combination A m A( mn ) A n is required. Yokoyama, 
Suyama & Tanaka supplied an integral expression from which this could be computed. In 
Ref . [8] this was extended to an explicit evolution equation for Nr mn -\ . Equivalent expressions 
were later given by Mazumdar & Wang [71]. 

Eq. (4.9) expresses the evolution of T 1 m with A. It can be verified that its evolution 
with A* satisfies an analogous equation 

Vn^T^ = —T I n W n m . (4.39) 

Here, the covariant derivative applies to tangent-space indices at A* and therefore operates 
only on m. The index / labels a tangent-space basis at A and is inert. Using (4.39) together 
with (4.31)-(4.32) we conclude 

T> m N n = -N m w m n: (4.40) 
^N f N {mn) = -N( rn) w T m - N( mr) w r n - N r w r {mn) . (4.41) 

Like Eqs. (4.9)-(4.10) these are covariantized versions of the evolution equations in flat field- 
space, using the correct coefficients w m n and w m [ r s) which appear in the Jacobi equation. 

Eqs. (4.40)-(4.41) should be solved by fixing A to be the late time of interest. During 
inflation the initial conditions would correspond to (4.25)-(4.26). One should then integrate 
backwards until A* corresponds to the horizon-crossing time for the fiducial scale fc*. 



5 Non-minimally coupled models 



In a multiple-field model with potential V and non-minimal coupling to the Ricci scalar, the 
action in the Jordan frame can be written as 



S 



d 4 x 



(5.1) 



where f((f) is positive definite but otherwise arbitrary. We distinguish Jordan frame quan- 
tities with a circumflex. The kinetic energy is X, assumed to be an arbitrary second-order 
combination of the field derivatives in the form 



(5.2) 



where Gij(<p) is the Jordan frame field-space metric. It has been shown that this action can 
be rewritten in the Einstein frame after a conformal transformation [25] 



The Einstein frame action is 

S 



I 



dV 



9fJLV fdllV- 



1 



where V = V / f 2 is the potential and 



1 



X = -■zGusTWdytfr' 



-M£R + X-V 



Gu = jG u + 



2 f,lf,J 

2~-pi p ' 



(5.3) 



(5.4) 



(5.5) 
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Cosmological observables themselves are not altered by this procedure [72-75]. However, one 
must be careful to restrict attention to clearly defined physical quantities; in particular, the 
curvature perturbation £ is not itself an observable [76]. 

The Einstein-frame action (5.4) is of the same form as Eq. (2.1). Therefore the results 
of §§3-4 are also applicable to models with non-minimal coupling. 

6 Conclusions 

In this paper we have computed the covariant three-point function near horizon-crossing for 
a collection of slowly-rolling scalar fields with a nontrivial field-space metric. After making 
a conformal transformation, this framework is sufficiently general to include scenarios where 
one or more fields are nonmimimally coupled to the Ricci scalar. The subsequent superhorizon 
evolution can be expressed using a version of the 'separate universe' approach. 

We concentrate on the broad class of models described by a <7-model Lagrangian C = 
X + V, where 2X = —Gjjd^ip 1 d^<p J , and obtain expressions for the two- and three-point 
functions. The presence of a nontrivial field-space metric leads to technical subtleties. First, 
to obtain a covariant formalism we must be careful to define perturbations as tangent-space 
vectors Q using the method of Gong Sz Tanaka. When computing correlation functions in- 
volving the Q , the Feynman diagram expansion introduces explicit factors of the 'trajectory 
propagator' which implements parallel transport along the inflationary trajectory. 

Second, new interaction vertices appear which involve explicit factors of the Riemann 
curvature tensor. Therefore the two- and three-point functions receive modifications of two 
types. The first follow from promotion of the flat field-space perturbation bp 1 to Q 1 and 
covariantize the result for G IJ = 5 IJ . As in general relativity, covariantization is achieved 
by exchanging partial derivatives for covariant derivatives and contracting all indices with 
the field-space metric. The second involve explicit factors of the field-space Riemann curva- 
ture. These modify both quantum interference effects operating near horizon exit, and the 
interactions between growing modes far outside the horizon. 

In §4 we developed a covariant version of the 'separate universe' formalism to account 
for these superhorizon-scale interactions. The covariant Jacobi equation (4.5) automatically 
incorporates curvature contributions which influence the evolution of the two- and three-point 
function. We have shown that this correctly reproduces the time-dependent growing modes 
near horizon-crossing generated by the apparatus of quantum field theory. In particular, it 
matches the two lowest-order divergences at single-logarithm order and the leading divergence 
at double- logarithm order. 

The Jacobi approach leads to covariant 'time evolution operators' T 1 m and T 1 ( mn ) which 
can be obtained straightforwardly by direct integration of (4.9)-(4.10). Together with the 
covariant gauge transformations derived in §4.3 these yield covariant l 5N coefficients' N m 
and iV(- mn ) which define the 'separate universe' expansion of the curvature perturbation £ in 
Eq. (4.30). This provides a clear and economical framework enabling perturbations to be 
evolved in a slow-roll inflationary model with nontrivial field-space metric. 

We always retain the option to abandon manifest covariance and work with the coor- 
dinate variation Sip 1 . The traditional separate- universe expansion for dip 1 is unchanged by 
the presence of a nontrivial field-space metric, and our predictions for the autocorrelation 
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functions of £ cannot vary because £ is a field-space scalar. The advantage of the covariant 
formulation is one of convenience and practicality. 

The phenomenology of the new Riemann-tensor terms may be interesting. In a canonical 
scenario, interactions among superhorizon modes are suppressed by three powers of 4>/H, 
and are therefore relatively slow. (Note that one should not regard this suppression as an 
indication of how large the bispectrum can become. Rather, it is an indication of the timescale 
over which it can evolve.) However, Eq. (3.17) shows that curvature-mediated interactions 
are suppressed by only a single power of <f)/H. In a model where the field-space curvature 
is 0(1) these could lead to much more rapid evolution. It will be interesting to study these 
effects in more detail, and we hope to return to this in future work. 
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A n-point functions in the in— in formalism 

The general theory necessary to compute n-point functions was set out in the papers by 
Maldacena [33] and Weinberg [46], and has been reviewed elsewhere [4, 61, 62]. In curved 
field-space this is complicated by the necessity to ensure that all quantities satisfy the correct 
tensor transformation laws. This means that the structure of the two- and higher n-point 
functions becomes more elaborate. Special factors, given by line integrals of the connection 
over field space, are required to ensure their indices reside in the proper tangent space. 

In this Appendix we briefly review necessary elements from the general theory, focusing 
on the changes necessary due to field-space curvature. 

A.l Two-point function 

After transforming to conformal time rj, defined by dt = a drj, the in-in generating functional 
can be written to quadratic order, 

Z = J [dQ ! + dQi] exp ^-i j d 3 xdr] a 2 Q/ ( A _ A ) Q J + <*- fn terms^ , (A.l) 

where = (QL,QL) and an overbar denotes matrix transposition. The time tq should be 
set well before horizon crossing of the modes under discussion, and r is the time at which we 



- 22 - 



wish to compute each correlation function. The ^-function terms have support at to and r, 
and enforce boundary conditions, to be described below. Finally, the differential operator A 
satisfies 

A/j = G u (v 2 v + 2^V V - did^j + a 2 M u , (A.2) 

where Mu is the mass matrix (3.11). 

We write the time-ordered two-point function between contravariant components of Q+ 
as .D++, 

D J + f\ V ,x;a,y) = (TQ^x^f (<7,y)>, (A.3) 

with analogous definitions for D |_, D_| and D In this appendix, we use unprimed 

indices to label the tangent space at time 77 and primed indices to label the tangent space 
at time a. Each two-point function is a field-space bitensor and a spacetime biscalar. The 
operator T denotes time ordering, with the convention that all ' — ' fields are taken later than 
all '+' fields and that time ordering on the ' — ' contour is in the reverse sense to the '+' 
contour. 

The rules of Gaussian integration enable us to calculate the D±±. They are obtained 
by inverting the quadratic structure A/j, 

K^UKif (-> 

Since an expectation value (O) inherits the tensor transformation properties of O, the asso- 
ciated two-point functions between covariant or mixed components of Q± can be obtained 
by raising or lowering the J and K' indices. To do so, one should use the metric evaluated 
at r\ or a, respectively. In (A.4) we have suppressed coordinate labels, but the differential 
operator acts only on rj and x. 

Mass matrix. For the remainder of this section, we ignore the mass matrix Aiu and treat 
each mode as massless. Small masses can be accommodated perturbatively if desired, but 
we will not do so in this paper. Note that this does not imply that we ignore all couplings 
between modes after horizon crossing: these are certainly important, because they describe 
how power is transferred from isocurvature perturbations to the adiabatic mode. These 
couplings will be retained when we discuss time evolution in §4. We are ignoring them only 
for a brief period around horizon exit. 

Tensor structure. First, consider -D++. The x and y dependence can be diagonalized by 
passing to Fourier space, 

D J 4 = J^Di«'(k)^<*-y\ (A.5) 
Neglecting the mass matrix as explained above, the mode function Di^_ (k) satisfies 

G I3 (v* + 2 a -V^ + k^j DiX = ~GfS(n - a). (A.6) 
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To factorize the tensor structure we introduce a bitensor II 1 k 1 which is required to solve the 
equation V^Ii 1 k 1 = 0, 

j _ dnV j &4> M N 

r > K ' ~ dry MN ~d~q~ ~ ^ ' 

The solution can formally be written as an ordered exponential, 

nV = Pexp (-^dr r 7 VAT»^-)Gf„ (A.8) 

where the integral is computed along the inflationary trajectory and the symbol 'P' denotes 
path ordering along it. Eq. (A.8) is the trajectory propagator introduced in §§3.2-3.3. It 
is simply the parallel propagator evaluated on the inflationary trajectory. In (A.8) we have 
chosen boundary conditions so that when 77 — > a the trajectory propagator satisfies H 1 k' — > 

The trajectory propagator allows the index structure in (A. 6) to be factorized. Taking 
= n Jft -'A + _|_ (where indices on Tl 1 k' are raised and lowered using the usual rules for a 
bitensor), it follows that the scalar factor A + _|_ satisfies the same equation as the propagator 
in flat field-space [1], 

Vl + 2-V v + k 2 ) A++ = -±-8(n-a). (A.9) 
1 a J a z 

The same factorization can be made for each Green 's function, so D££ = U JK ' A ±± . With 
vacuum boundary conditions, the 5-function terms at ro in (A.l) require A + _|_ to be approx- 
imately positive frequency there, and A |_ to be approximately negative frequency. The 

^-function terms at r require A ++ (r, a) = A y(r,o~) f° r an 17 ■ Finally, A and A^ are 

the Hermitian conjugates of A ++ and A |_, respectively. 

If the initial conditions at ro correspond to the vacuum, then D ++ should be approxi- 
mately positive frequency at that time. If ro is well before horizon-exit, then 

i-r^l /1 ^J'/, w , n \3r/, , \ttIJ' h * ( (I - ikri)(l + ika)^^'^ ri < a 
(TQi( kl , v )Qi (k 2 , a)) c (2vr^( kl + k 2 )Il" ^ x | ^ + ^ _ ^ , 

(A.10) 

where k is the common value of |ki| and | ki2 1 , and this estimate is valid for values of rj and a 
within a few e-folds of horizon-crossing at \krj\ = \ka\ = 1. A subscript '*' denotes evaluation 
precisely at the horizon-crossing time. In Eq. (A. 10) the trajectory propagator plays an 
essential role in ensuring that the right-hand side has the correct bitensorial transformation 
law. 

Now consider D |_. Imposing the boundary condition that D |_ is approximately neg- 
ative frequency at the initial time, and equals D ++ at time r, we find 

(TQLQli,v)Q+Ql2,<t)) = (2r) 3 5(k 1 + k 2 )n"'||(l + ik V )(l - ika)e ik ^\ (A.ll) 
A. 2 Three-point function 

Each term in the third-order action (3.14) makes a contribution to the three-point correlation 
function (Q^Q^ Q^ 3 )j or equivalently . Vertices are constructed from two copies of the 
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action, one for each of the '+' fields and ' — ' fields of Eq. (A.l). Vertices for '+' fields appear 
with a factor of +i; vertices for '— ' fields appear with a factor of — i. We apply Wick's 
theorem to produce all ways of pairing indices, using G ++ of (A. 10) to pair two '+' indices 
and its complex conjugate G to pair two '— ' indices. We use G h in (A. 11) or its complex 
conjugate to pair a mix of '+' and ' — ' indices. Finally, we integrate over all possible 
spacetime positions for each vertex. 

Trajectory propagator. In curved field-space, this procedure is modified by the appear- 
ance of the trajectory propagator in each two-point function. Consider a typical interaction 
appearing the the third-order action, such as the term 



/2 
d 3 xdr? { - —^—^ < j ) I 'Q II V v Q J 'v v Qj,y (A.12) 



To proceed we should expand background quantities around a reference scale k*, in the 
direction of the inflationary trajectory. This can be determined by analogy with the Taylor 
expansion (2.5), which was constructed along a geodesic. To separate the various tangent- 
space indices which appear, we adopt the following conventions: tangent-space indices at the 
time of observation, r, are labelled I, J, K. (In the main text, the time of observation is 
expressed as an e-folding number N .) Tangent-space indices at the time of horizon-crossing 
for the reference scale k* are labelled i, j, k. (In the main text, the time of horizon-crossing is 
expressed as A 7 *.) Finally, indices associated with the integration variable rj are given primed 
indices I', J', K'. The background factor <f)/H gives 

n I ' i (^+V N ^N + ---) , (A.13) 



H{ V ) \H ' " H 

where A" = — In ] fc* 77 ] represents the number of e- folds since horizon exit of the reference 
scale. We should make an analogous expansion for the metric which is used to contract the 
two copies of T> V Q J ' , 

G J ' K ' = n J ' j U K ' k Gi k . (A.14) 

In this case, higher-order terms in the Taylor expansion vanish because the metric is co- 
variantly constant. A final example is the Riemann tensor, whose series expansion will be 
required to compute the higher-order contribution (A. 22), 

R i>(j'K')L>h. = u I ' i IL J ' j IL K ' k ( + V N R^ l ^N + ■■■) . (A.15) 
H V H H J # 

One can see that the three factors of the trajectory propagator in Eq. (A.15) will be 
common to all terms in the three-point function. These propagators carry dependence on 
the integration variable r\ as well as the reference time A*. However, each two-point function 
which connects an external field with a field at the vertex will introduce a propagator factor 
H 1 1>, which is a function of r\ and N. Therefore all 77 dependence cancels, leaving a propagator 
that relates the two tangent spaces at A^ and A^. This may be factored out of the vertex 
integral. The remaining details of the calculation correspond to those in flat field-space. 
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Results. Three of the contributing terms are covariantized versions of those present in flat 
field-space: 



(I + |) + cyclic, 



2M£ H* 



(A.16) 



VQidiQ J diQj 



1 iff 



G{ K (k 2 -k 3 ) 



k 2 k\k 2 k 3 

h~ k 2 



+ cyclic, 



(A.17) 



2AfJ-ET 

A ijk D 



* d k 



2M 2 ,H, * 



1 k 



1 , ^3 



+ cyclic. (A. 18) 



These expressions are valid under the same conditions as the power spectrum (3.12): we 
must wait sufficiently many e-folds that decaying power-law terms have become negligible, 
but Eqs. (A.16)— (A. 18) are unreliable when N = — In | fc^-rj 3> 1. These three contributions 
can be combined to give 



1 <P* r jk 



2 7.2 



/Co Jv 



2^3 



+ -h{k 2 -k 3 ) 



+ cyclic. 



(A.19) 



We now consider the remaining two terms involving the Riemann tensor: 
2a 3 



-R IiJK)L 0LV v QiQjQ K 



4 

— j 

3 



j^ijk -j ** j^i(jk)m<Pm 



AT i 1 M I. 1.2, ^1^3 

fef ( 7e - A + In — j - k t k x H — — 



+ cyclic, 



(A.20) 



• ^R {I]LMIJ -' K) ^mQiQjQk 
6 



3 



kt 



1 



fc? A - In — - 7E - - + -kf - ktK 



9 



+ cyclic. (A.21) 



Finally we consider the next-order correction to the K l{ ^ JK ^ L ^ term: 



— R HJK)LlM^ MNVriQlQjQK 



^ijk ^ j^i(jk)m;n ( Pm Qn 

3 H* 



kfN 2 + kf[ 1 2 --+ln^ 

-2A; t A; 1 2 (ln^+7 E -l) + 5 
+ cyclic. 



7T 



(27E 



+ ln 



k 2 k 2 k 3 r kt 



In — + 7e 



(A.22) 
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The presence of terms proportional to positive powers of N, which spoil the validity of these 
formulae when N S> 1, is explicit in Eqs. (A.20)-(A.22). 
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